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NOTATION 


One-half parallel length of model 
Maximum beam of model 
One-half maximum beam of model 


Residual resistance coefficient 


Calculated wave resistance coefficient 


Froude number 
Depth of immersion 
Acceleration due to gravity 


Draft 


Struve’s function 
Total length of model 


Parabolic length 
f 1 sin yéd é 
0) 
fp ecos vas 
Ratio of parallel length to total length 


Resistance as calculated from Michell’s integral 
Wave resistance coefficient 

Wetted surface 

Velocity 

Bessel function of the second kind 

Surface elevation due to symmetrical disturbance 
Surface elevation due to wave system 


Waterline coefficient 


Displacement 
Block coefficient 


Nondimensional vertical measurement 


Qq > Mm > 


il 


Wave length 


Nondimensional longitudinal measurement 


Density 


Prismatic coefficient 


INVESTIGATION OF WAVE EFFECTS PRODUCED BY A THIN BODY - 
TMB MODEL 4125 


by 


Georg P. Weinblum, Janet J. Kendrick, and M. Allison Todd 


INTRODUCTION 


An opportunity to check experimentally Michell’s wave-resistance theory! arose when 
a friction plane or ‘‘thin body’’ was constructed for investigations of the frictional resistance 
of painted surfaces. In his wave-resistance theory, the only assumption made by Michell as 
to the form of the ship was that the inclination of the tangent plane at any point of its surface 
to the vertical median plane should be small. 

Although the TMB model is quite narrow, it cannot be considered as a thin plank; thus, 
appreciable wave effects must be expected when the model is towed at medium and high 
Froude numbers. 

- The dimensions of the body are shown in Figure 1. It has a parallel middle body and 
parabolic ends with vertical sides. Apart from the rectangular form of the sections, this body 
forms an ideal Michell’s ship as defined above due to its very low beam-length ratio 
[B/L = 0.0265] and beam-draft ratio [B/H = 0.183]. 


Figure 1 - Shape and Dimensions of Friction Body, TMB Model 4125 


L = total length = 252 inches H = draft = 36 inches 
66 inches a= half parallel length = 60 inches 
0. 4762 


Ul = parabolic length 


B=maximum beam = 6.69 inches p= 2a/L 
b = half beam 3.345 inches 


IReferences are listed on page 14. 


It was decided to calculate the wave resistance experienced by the model and the 
wave profiles along the hull and to compare these theoretical values with the results of ex- 


periments. 
Auxiliary resistance integrals, which simplify the work considerably, were not avail- 


able at the time. Since, however, the draft-length ratio H/L is not covered by the systematic 
computations now being made at the Taylor Model Basin, an independent evaluation was need- 


ed in any case and justified the appreciable work involved. 


RESISTANCE 
The original form of Michell’s resistance integral may be reduced to the dimensionless 
form? 
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15: 1, represent the symmetrical and antisymmetrical parts of the hull with respect 
to the midship section, and 


” = speed of advance. 


Since for a symmetrical hull /*(y) = 0, Equation [1] reduces in this case to 


Goods [4] 


The intermediate function J*(y) is given in this case by 
I*(y)= —2E,|sinypM, {7 (-p)} + cos ypM, {y(1-p)} | [5] 
as can be seen by substituting the form of the waterline equation 7, in Equation [3]. 
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Graphs of E? over 4 and of sin yp M{ {y(- p)\ + cosy pM, {yap} over y were plot- 


ted; from these the necessary values could be obtained. 
The integral was evaluated at the singularity y/y)= 1 by integrating over a narrow 


range between + y, (1 + €) where « < 0.01 using the formula 
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was evaluated from the graphs of the integrand by means of a planimeter. The remainder 
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was calculated from expansion of the integrand. 


J’ (y)dy 


The resistance was then calculated for two additional drafts, H = 20 inches and 10 
inches, and at two Froude numbers for an infinite draft. 
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Havelock has given an example showing how to correct resistance data obtained from 
friction body experiments for the wave resistance.* In the discussion of this paper, Weinblum 
has pointed out that these estimates may be much too low for moderate Froude numbers since 
they are based on a parabolic waterline while the actual form can be much fuller (see page 
269, Reference 3). 

The results of the computations are presented in various ways. Figure 2, reproduced 


from Reference 4, shows the comparison of the calculated wave resistance coefficient 


= ug with the ‘‘residual resistance’’ coefficient C.. Schoenherr’s frictional 


+ Su? 
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resistance line was used to evaluate Che 


The coefficients reach their maximum at moderate Froude numbers which is character- 


istic of very full forms. Figure 3 shows the calculated wave resistance coefficients 
T= aia Tek plotted against y= | for three depth-length ratios. 
8p9 B?H? QF? 
eas 
Figure 4 shows the wave resistance coefficients of a parabolic form (for which p = 0) 
with the same proportions as the friction body. It can be seen that for this form, the coeffi- 
cient of resistance is much lower than that of the fuller form except in the region of the large 


hump. 
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Figure 3 - Wave Resistance Coefficients r Plotted Against Speed 
for Three Different Drafts 
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Figure 4 - Comparison of Wave Resistance Coefficients 
of Parabolic and Full Forms 


To obtain a better basis of comparison between theory and experiment, it was suggest- 
ed that the influence of the sharp corners on the resistance should be eliminated by towing 
the model at a shallower draft H = 20 inches. Then the measured and calculated differences 
in resistance for the two drafts could have been calculated. However, it was not feasible to 


tow the model at a shallower draft at the time the tests were run. 


WAVE PROFILES 


The work done on the calculation of the wave profiles of the friction body was based 
upon information given by W.CS. Wigley.° His investigations were performed upon a body 
having a short parallel middle body, 1/16 of the total length, and parabolic ends. 

Since the subject of wave profiles has not been treated as frequently as that of wave 
resistance, the computations are presented here in detail. 

Two formulas given by Wigley in dimensional form were used in the calculations of the 
wave profiles. The first gave the surface elevation due to wave making and the second the 
symmetrical disturbance. 


In dimensionless coefficients these become 
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where Z, is the surface elevation due to non-wave portion of disturbance caused by motion of 
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the form,* and Z_ is the surface elevation due to wave system. 
The nondimensional longitudinal distance €is positive when measured in the astern 


direction. The basic functions P and @Q are defined as follows:7’ ® 
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In the present calculations, the values of the P functions were obtained from a graph prepared 


by Professor Lunde.” 
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where H,(¢) is Struve’s function and Y,(¢)is the Bessel function of the second kind. The 
values of Q)(u) may be calculated for values of ¢ < 16 using the tables of H,(¢) and Y(t) 
given by Watson.° For values of ¢> 16 (for which the values of H(t) and Y,(t) are not given), 
the following method was used: 


When ¢ is large (in this case for ¢ > 16) 
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*<¢The function Q does not oscillate but is monotonic, and the terms in Q represent,a symmetrical disturbance 
of the surface in the neighborhood of the form which dies away quickly both fore and aft of the form. It absorbs 
no energy owing to its symmetry and therefore does not affect the resistance.’ 6 
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TABLE 2 
Values for Q,(w) ={- Qo(t) dt 


Q,(u) | u @,(u) 
24.988 | 20..5 68.636 
26.775 | 21.0 70.803 
28.577 | Zio 8 72.963 
30.411 | 22.0 75.153 
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Evaluating to the first two terms only 
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Q,(u) was evaluated for 16 <¢ < 30 using Equation [8]. A tabulation of these values is given 
in Table 1. 

The values of Q,(u) = i Q(t) dt were obtained from these by integration, see Table 
2. 

The P-functions are not defined for negative values of ¢. P(—|¢|) = 0 was used. Thus, 
it can be seen from Equation [6] for Z, that for 0 < € [1 —p only the terms P,(y€) and 
Rev? (y €) are defined. For 1 -—p = €S1 +> the term Peet ty(€-1-=p)}is defined. For 
1+ pSé&S2 the term Pra {y(€ — 1+ p)} is also defined and for €22 all are defined. It 
would appear then that the first two terms mentioned would give the ordinates of the bow wave 
system, the next the ordinates for the waves proceeding from the foreshoulder, the fourth the 
ordinates of the waves proceeding from the aftershoulder, and the last two the ordinates of 
the stern wave system. 

In the example given in the present paper (Figure 6), these four waves are plotted 
separately to show the components of the total wave system contributed by each of the terms 
considered above. 
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The Q-functions appearing in the equation for the symmetrical disturbance’ are defined 


for negative values of ¢ as: 
Q,(-t)=Q,(t) and Q,(—-t) = —Q,(t) 


Thus all the terms of Equation [7] are defined for all values of €. 
The evaluation of Equations [6] and [7] give the wave profiles in terms of dimension- 


less coefficients containing Z, and Z). The multiplication of these values by the constants 
ies 2b 
m™y(1—p) m™y(1—?p) 

Figure 7 the calculated and observed profiles are plotted to an inch scale. in Figure 6, how- 


an respectively, transforms them into the dimensions of 6. In 


ever, the components are plotted in dimensionless form. 

The comparison was made as outlined above for five different speeds. At the lower 
speeds (4.5 and 5 knots) the agreement between the observed and calculated profiles is good; 
for higher speeds the agreement of the crests is good but the troughs of the calculated pro- 
files are deeper than those observed (Figure 7). This is at least partly accounted for by the 


assumption of infinite draft made throughout the calculation. The influence of a finite draft 
will be more marked at higher speeds. This will be evident from a consideration of the 
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Figure 5 - Wave Resistance Coefficient C, Plotted Against Froude Number 
for Parabolic and Full Forms 
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Figure 6a - Main Surface Elevation 
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Figure 6b - Symmetrical Disturbance 


Figure 6 - Wave Profiles for Speed of Advance of 6 Knots Plotted in Dimensionless 
Coefficients and Showing Components 


12 


V=4.5 knots 
| | 


ts Be 
SS Se -4 
=0.326 V=5 knots 
yee -8 


we 


= 
7 


- 
- 
- 


V=7.72 knots 


e—e Experimental Wave Profile ee co) 
+----+ Calculated Wave Profile ee 
a Additional Observed Profiles abi laolweecea © Si aaa 4 
Destine aca 
Sela loka rae 
-8 
ai I) 
-12 
1.0 (o) -1.0 
Stations 
FP 
a | We | ye | 
20 15 10 5 (0) 
Scale for Observed Profiles 


Figure 7 - Comparison of Calculate 


d and Measured Wave Profiles 


for Five Different Speeds of Advance 
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exponential term relating the below-surface disturbing pressures to the Ee length at the sur- 
a2mt 


A 


face and the depth of immersion. These forces vary approximately as e where fis the 


depth of immersion and ) is the wave length. Thus as A increases with increasing speed, 
— 2af 


e % also increases. The effect upon the water surface elevation of the finite draft of the 
model will, therefore, be greater at higher speeds than at low, and better agreement than that 
actually obtained cannot be expected. 
It was found that measuring the wave profile at the bow sometimes presents serious 
difficulties because of the steep slope of the wave in a transverse plane. Photographs of 


wave contours along the ship may therefore prove to be rather unreliable. 


CONCLUSION 


The qualitative and quantitative agreements between the calculated wave resistance 
coefficient C, and the residual resistance coefficient C,(as shown in Figure 2) are very 
satisfactory. In particular, prediction of the ‘‘third’’ hump* was possible in the resistance 
curve of Figure 3 at F = 0.24 which at first escaped the attention of the experimenters. 

Because of this good agreement, the following deductions from the computations can 
be relied upon without further checks: 

It is well known that because of the linear character of Michell’s analysis the relation 
R ~ B? holds. 

It can also be demonstrated that at very large Froude numbers (y > 0) the integral r 
becomes independent of H/L provided this ratio is finite. Thus the three curves on Figure 3 
must finally coincide and R ~ BH, i.e., the wave resistance in this range is proportional 
also to the square of the depth H. However, it can be seen from the figures that this asymp- 
totic relation does not hold for high-speed displacement ships. For example, at a destroyer 
speed of F = 0.6, y= 1.5. In this case R is proportional to H1:5 only; at F = 0.25, R~H? 
where 7 is already less than unity. 

The studies of the method of formation of wave systems described in this report ex- 
plain to some extent the characteristic form of the wave-resistance curve. 

It can readily be seen that with increasing speed and therefore with increasing wave 
length, the crests of the bow wave system will alternately reinforce and dampen those of the 
stern and shoulder wave systems. Similarly, combinations of these systems will change with 


changing speed. These reinforcing and damping effects of the wave systems account for the 


*The hump in the resistance coefficient curve occurring at F = 0.5 is here denoted as the ‘‘first’’ hump, at 
F =0.3 as the ‘‘second’’ hump, etc. (see Figure 3), contrary to the custom in naval architecture by which the 
hump at the highest speed is called the last hump. This change appeared to be necessary since from a mathemat- 


ical viewpoint there are an infinite number of humps between 0.5 2 F 2 0 (see also Reference 2, page 23). 
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peaks and hollows seen in the curve of wave resistance plotted against speed. However, 
after a certain speed is reached, even the second wave of the bow wave system will be astern 
of the ship, and there can be no further reinforcement of the other systems. Above this speed, 
therefore, the wave resistance curve will have no further peaks but will show a steady fall 


with further increase in speed. 
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